Josephson Plasma Resonance as a Structural Probe of Vortex Liquid 
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Recent developments of the Josephson plasma resonance and transport c-axis measurements in 
layered high T c superconductors allow to probe Josephson coupling in a wide range of the vortex 
phase diagram. We derive a relation between the field dependent Josephson coupling energy and 
the density correlation function of the vortex liquid. This relation provides a unique opportunity 
to extract the density correlation function of pancake vortices from the dependence of the plasma 
resonance on the afe-component of the magnetic field at a fixed c-axis component. 



Thermal suppression of Josephson interlayer coupling 
by vortex fluctuations pi is the key factor, which en- 
hances vortex mobility and limits useful applications of 
Bi- and Tl-based high T c superconductors (HTSC). Re- 
cent discovery of the theoretically predicted J|,H Joseph- 
son plasma resonance (JPR) effect in the mixed state of 
highly anisotropic layered high temperature and organic 
superconductors [|]-^| provides a new and unique oppor- 
tunity to probe the Josephson coupling of layers in a wide 
range of magnetic fields and temperatures. An alterna- 
tive way to probe Josephson coupling in the vortex state 
is direct c-axis transport measurements [ flO[ . 

A minimum necessary condition for the existence of a 
narrow JPR line is short-time/short-range phase coher- 
ence between the superconducting layers. In the Lon- 
don regime this coherence is quantitatively character- 
ized by the spatial and temporal behavior of the "local 
coherence parameter" C n (v,t) — cos[y ni „+i(r, t)], with 
ip n ,n+i = <Pn+i ~ <Pn - (2tts/^ q )A z being the gauge- 
invariant phase difference between layers n and n + 1. 
Here r = (x,y) and z are coordinates in the ab plane and 
along the c axis, s is the interlayer spacing. 

The first important condition is that the typical time 
of variation of this factor (phase slip time) must be much 
larger than the period of plasma oscillations. Their ratio 
was estimated |nj as ~ 10 4 at T — 45 K and magnetic 
field B z = 0.1 T. In this case the plasma mode probes a 
snapshot of the instantaneous phase distribution. 

In the decoupled liquid state the cosine factor C n (r, t = 
0) also rapidly oscillates in space because correlations 
between pancake positions in neighboring layers are al- 
most absent. The possibility of a sharp resonance in the 
situation when fluctuations of local Josephson coupling 
are much stronger than its average is an interesting and 
nontrivial issue. JPR resonance in such situation occurs 
because small phase oscillations induced by the external 
electric field change slowly in space and average out these 
rapid variations [ p~2]JT^ ] . 

In this Letter, using the high temperature expansion 
with respect to the Josephson coupling, we establish a re- 
lation connecting the plasma frequency with the density 



correlation function of the pancake liquid. A similar re- 
lation was derived in Ref. |Q assuming the Gaussian 
distribution for the interlayer phase fluctuations. We 
demonstrate that the pair distribution function of the 
liquid is unambiguously connected with the dependence 
of the plasma frequency on the in-plane component of 
the magnetic field for fixed c-axis component. Therefore, 
this dependence can be used to extract quantitative in- 
formation about the structure of the liquid phase. 

The essential physics of JPR is captured by a simpli- 
fied equation for small oscillations of the phase difference 
ip' n n+ i(r,uj) induced by the external microwave electric 
field with the amplitude V z and the frequency uj (see 
Ref. @) 
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where lu (T) — c/ v / e^A c (T) is the zero field plasma fre- 
quency, e c is the high frequency dielectric constant, A c 
and X a b are the components of the London penetration 
depth, Ej — Eq/Xj is the Josephson energy per unit 
area, Eq = ,s$Q/167r 3 A^ b , Xj = js is the Josephson 
length, and 7 = X c /X a b is the anisotropy ratio. The 
inductive matrix L is defined as LA n — ^ m L nm A m 
with L„ m s» (A q6 /2s) exp (-\n - m\s/X ab ). We neglect 
in Eq. (Q) time variations of C„(r,t) assuming them to 
be small during the time l/u>. Eq. (Q) also does not take 
into account charging effects [H and quasiparticle dissi- 
pation. These effects don't influence much the position 
of the resonance but may modify its lineshape. 

Static configurations (p„ tn +i(r) are mainly determined 
by thermal fluctuations of pancake vortices and quanti- 
tatively characterized by the average value of the cosine 
factor, C(B,T), 



C(B,T) = (cos^», n+1 (r)) 
and its static correlation function 5(r, B,T), 
5(r,B,T) = (cos[<£ n , n+ i(r) - y>n,n+i(0)]), 



(2) 



(3) 



where (. . .) denotes the thermal average. In this Letter 
we consider the pancake liquid state where the correla- 
tions between pancake arrangements in different layers 
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are almost absent. In such situation C„(r) rapidly oscil- 
lates in space, so that S(r, B,T) drops at distances of 
the order of the intervortex spacing a = (Qq/Bz) 1 / 2 , and 
C(B,T)<1. 

An important observation is that in spite of rapid vari- 
ations of C„(r) the oscillating phase tp' = ip' n „ +1 (r) varies 
smoothly in space. The typical length scale L v of <p' 
variations can be estimated by balancing the typical ki- 
netic energy of supercurrents, E (ip') 2 / 1? , with the typ- 
ical Josephson energy, Ej(ip') 2 a/ L v , in the vortex state 
which is strongly disordered along the c axis. This gives 
L v = \ 2 /a. Smoothly varying tp' effectively averages 
rapid variation of C n (r) and the plasma frequency is sim- 
ply determined by C(B,T) (see Refs. 



<p n , n+l (r) = ¥>S + i( r ) + <n+i( r ) " (2vr S /$ )S,y. (9) 

The phase <p^ n+1 {Y;Y nv ,r n+ i tU ) is the singular part of 
the phase difference induced by vortices at positions r n „ , 
rn+i,i/ when Josephson coupling is absent [Ej = 0): 



u; 2 p (B,T)=u; 2 (T)C(B,T)/C(T), 



(4) 



Here the factor C(T) = C(0,T) takes into account the 
suppression of zero field plasma frequency by phase fluc- 
tuations (see below) . Fluctuations of C n (r) smoothened 
over the area 1? , C{L V ) = L~ 2 J r<L drC n (r), produce 
inhomogeneous broadening of the JPR line. Calculating 
the mean squared fluctuation of C{L V ), ([C(L V ) — C] 2 ) w 
a 2 /L 2 = a 4 /X 4 j, we obtain the estimate for the inhomo- 
geneous line width 



5(u, 2 )« W 2 (7> 2 /A 2 . 



(5) 



Detailed calculations of the intrinsic lineshape due to this 
mechanism will be published elsewhere p3[ . 

We establish now a relation between plasma frequency 
and the vortex liquid structure. In general, the Josephson 
coupled system in the London regime, B z <C H C 2, may 
be described in terms of vortex coordinates r„„ (index v 
labels vortices in the layer n) and by the regular (spin- 
wave type) phase difference ip^ n+1 (r). The free energy 
functional in terms of these variables is p2| 

Here T v (y nv ) accounts for the two-dimensional energy 
of pancakes and also their electromagnetic interaction in 
different layers jl6|]. The second term, 
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is the energy of intralayer currents associated with regu- 
lar phase fluctuations, and the third term, 
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is the Josephson energy. We split the phase difference 
tp n ,n+l into vortex part, regular spin-wave part, and con- 
tribution coming from the in-plane magnetic field B x , 
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where 4> v {v) is the polar angle of the point r. 

We calculate C(B, T) with the functional (^) using high 
temperature expansion jll[ with respect to Tj as: 



C(B,T) « §/ drS(r,B,T) = /(B, T) 



EnB 



2TB 7 



(11) 



a 2 J drSJr, B, T) is a universal function 



Here/(B,T) 

and Bj = <f> /\ 2 j. Eq. @ is valid until C(B,T) < 1 
which corresponds to the field and temperature range 
TB Z > E Bj. Comparing Eqs. @ and Q with 
Eq. (^) we obtain an estimate for the relative line width 
in the liquid state due to the inhomogeneous broadening, 
duip/ujp w T/Eo <C 1. 

In the lowest order in Ej spin-wave and vortex degrees 
of freedom do not interact: 



S(r,B) = S v (r,B z )S r (r)cos(2nsB x y/<Z> ), 



(12) 



where the correlation function S v (r) is determined by the 
functional J- Vl while S r (r) is determined by T v . In the 
pancake liquid regime, which we consider in this Letter, 
B x penetrates almost freely into the sample and has no 

(v) (r) 

influence on the phase fluctuations fnn+i and p n ' n+ i- 
An extra spatial dependence of ip n , n +i induced by B x 
gives possibility to probe phase correlations at given B z 
by measuring i^-dependence of the plasma frequency. 

We now express the vortex phase correlation func- 
tion S v (r) via the density correlation function of pan- 
cake vortices induced by B z . The difference $ u (r) = 
</>„„ + i (r) — <^„„-i_i(0) can be connected with pancake 
densities /0„(R) = ^(-^ — Vnl/ ) m ^ ne l & y ers as 



*„(r) = j dR[ Pn (R) - Pn+1 (R)}/3(r,R) 



(13) 



Here /3(r, R) = <f> v (r/2 - R) - <f> v (-r/2 - R) is the angle 
at which the segment connecting points —r/2 and r/2 is 
seen from the point R, 

cos/?(r,R) = (i? 2 -r 2 /4)[(i? 2 +r 2 /4) 2 - (Rr) 2 ]- 1 / 2 . 

The function /3(r, R) has a jump of 2ir when point R in- 
tersects segment [—r/2, r/2]. Using the Gaussian approx- 
imation for (r) we obtain for vortex phase correlation 
function S„(r) = (cos[$„(r)]) 



S v (t) =exp[-F„(r)] 



(14) 
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with 



F v (r) = ([$„(r)] 2 )/2. 



(15) 



Using Eqs. (|13|) and ( |15[ ) we connect ^(r) with the pair 
distribution function h(r) defined by the relation 

n 2 v h(r) = (p„(0K(r)> - (p n (0)p„ +1 (r)) - n^(r) (16) 



Fv(r) = 

J dRdr'h(r') 



;(r,R+|J -/?(r,R-| 



(17) 
n 2 



For this identity n v J drh(r) = —1 has been used, n v = 
(p n (r)) = aT 2 . Performing integration with respect to R 
and averaging over directions of r' we finally obtain 

F v (r) = -mil J r'dr'h(r')J(r,r'), (18) 

where the universal function J{r, r') is given by 

J(r, r') = 2irr 2 [ln(r'/r) + 2] , r > r, (19) 
J(r, r') = 4tt [2rr' - r' 2 - (r' 2 /2) bx(r/r')] , r' < r. 

Eqs. (|lj), (|lj) and ( |l9| ) represent general relations con- 
necting phase fluctuations with vortex density fluctua- 
tions. The function F v (r) has the following asymptotics 



F v {r)^-8ir 2 



dxx 2 h(x), r ^> a, 



nr 2 (a [°° \ 

F v (r) « — I In - + 2 - 2tt y dxxh(x) Inxj 



r <C a. 



Linear increase of (r) at large r indicates exponential 
decay of phase correlations at large distances. 

The problem of calculating of S v (r) is now reduced to 
calculation of the integral dig ) with the density correla- 
tion function of the liquid. This function is not available 
analytically. It can be calculated from Monte Carlo sim- 
ulations or (approximately) using the density functional 
theory JTa]. In the decoupled pancake liquid regime the 
correlations between two dimensional pancake liquids in 
different layers are very weak because the energy of the 
magnetic interlayer interaction of disordered pancakes 
has an additional small parameter s 2 /X 2 b in compari- 
son with the intralayer energy. If these correlations are 
neglected then the vortex system is equivalent to the one- 
component two-dimensional Coulomb plasma which was 
studied extensively Jig] . The pair distribution function 
^2_o( r : B Z1 T, Eq) within this approximation has an ex- 
act scaling property /z^d^ B z ,T,E ) = /i2£>(?"/ a j T/Eq), 
i.e., it depends on magnetic field only through the spa- 
tial scale. As follows from Eq. ( O ) this scaling prop- 
erty is also transferred to the function F v (r), F v (r) = 
F v (r/a, T/Eq). If, in addition, fluctuations of the regular 



phase are neglected then at B x = the function /(T/Eq) 
in the expression ( ftl| ) for C(B,T) is field independent. 
Therefore, the scaling property C(B Z ,T) oc 1/B Z is a 
consequence of the approximation of completely decou- 
pled pancake liquids in pinning-frcc layers when spin- 
wave phase fluctuations are negligible. 

We generated the pair distribution functions h(x) with 
x = r/a for the two-dimensional pancake liquid using 
Langevin dynamics simulations. The model and the al- 
gorithm have been described in Ref. |L9| . Fig. 1 shows an 
example of the pair distribution function h{x) for tem- 
perature T = 0.0127ri?o (~ 1.7 T m ) and the correspond- 
ing function F v (x) obtained by numerical integration of 
Eq. (|l^) . To show a weak waving of F v (x) due to liquid 
correlations we also plot the derivative of F v (x). 

As temperature approaches T c the role of fluctuations 
of the regular phase progressively increases. The spin- 
wave phase correlations decay algebraically 



S r (r) = (U/ry a , a = T/2nE (T), 



(20) 



where £ a b(T) is the superconducting correlation length 
[1/Cah determines the upper cut-off for momenta char- 
acterizing spatial variations of <Pnli+i( r ) in Eq. (j7|)]. As 
follows from Eq. ( Jl2| ) the contribution from regular phase 
fluctuations becomes essential if S r (r ~ a) <C 1. This 
gives the condition T > E Q (T) ln(H c2 /B). 

To relate u> p (B,T) with the plasma frequency at zero 
field and the Josephson coupling Ej oc A~ 2 we have to 
take into account that these quantities are also renor- 
malized by the phase fluctuations. Their suppression is 
determined by the cosine factor at B — 0, C(T), which 
was estimated in Ref. iGOf as 



C(T) = (U/^jY 



(21) 



Finally using Eqs. (|), @, @), @, and @ we ob- 
tain an expression for plasma frequency which incorpo- 
rates effects of vortices and of regular phase fluctuations 
at all temperatures : 



co 2 (B,T) Eq(T) 
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(22) 
(23) 



where Jo(x) is the Bessel function. Note that regular 
phase fluctuations reduce the power index in the B z de- 
pendence of C(B Zl T) at B x = and make it temperature 
dependent. Using the typical parameters for optimally 
doped Bi 2 Sr 2 CaCu20 8 _ x [A o6 (0) = 2000 A, T c = 90 K] 
we obtain that the index a increases from a ~ 0.05 at 
T = 50 K to a w 0.18 at T =75 K. 

Measurements of independence of the resonance fre- 
quency uj p at fixed B z and T allow one to obtain, in prin- 
ciple, the pair distribution function h(x,T), i.e., to ex- 
tract quantitative information about correlations in the 
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vortex liquid. Using Eq. (|22|) the function f(b) can be 
found from the dependence lo p (B x , B z ) as: 



(24) 



With f s (b) known, the function F v (x) may be found by 
the reverse Fourier transform from Eq. (pj 



F v {x) = - In 



x Za / bdbf s (b)J (bx) 
o 



(25) 



Finally, the pair distribution function h(x) can be ex- 
pressed through F v (x) using Eqs. ( |l8|) and (|Tq|) : 



h(x) 



1 d 

8ttx 3 dx I efa; 



d 3 F v (j 
dx 3 



(26) 



Unfortunately, high order differentiations in the last 
equation complicate practical realization of the proce- 
dure. As it is shown in Fig. 1, oscillating behavior of 
liquid correlations is hardly seen in F v (x) but the first 
derivative, F^[x), may expose it quite clearly. The ex- 
perimental angular dependence of the plasma frequency 
HU was found to be in a good qualitative agreement 
with Eq. (p2|). However an accuracy of existing measure- 
ments is not sufficient to perform the described quanti- 
tative analyses. 

In conclusion, we obtain a general expression connect- 
ing the plasma frequency with the density correlation 
function of pancake liquid. We demonstrate that the de- 
pendence of the plasma frequency on the magnetic field 
component parallel to the layers B x at fixed perpendic- 
ular component B z contains full structural information 
about the vortex liquid at given B z and temperature. 
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FIG. 1. (a) Pair distribution function h(x) obtained by 
Langevin dynamics simulations of two dimensional vortex liq- 
uid for T/ttEo = 0.012; (b) Plot of mean squared phase dif- 
ference between points separated by distance r = xa, F v (x), 
[see Eq. (Jl5|) ] obtained by numerical integration of Eq. ( |lq ) 
with plotted above h(x). Derivative of F v (x) is also plotted 
to enhance weak oscillations due to liquid correlations. 
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